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ABSTRACT 

Context. Warm debris disks are a sub-sample of the large population of debris disks, and display excess emission in the mid-infrared. 
Around solar-type stars, very few objects (~ 2% of all debris disks) show emission features in mid-IR spectroscopic observations, that 
are attributed to small, warm silicate dust grains. The origin of this warm dust can possibly be explained either by a recent catastrophic 
collision between several bodies or by transport from an outer belt similar to the Kuiper belt in the Solar System. 
Aims. We present and analyse new far-IR Herschel/PAcs photometric observations, supplemented by new and archival ground-based 
data in the mid-IR (VLTI/Midi and VLT/Visir), for one of these rare systems: the 10-16 Myr old debris disk around HD 113766 A. 
We improve an existing model to account for these new observations. 

Methods. We implemented the contribution of an outer planetesimal belt in the Debra code, and successfully used it to model the SED 
as well as complementary observations, notably Midi data. We better constrain the spatial distribution of the dust and its composition. 
Results. We underline the limitations of SED modelling and the need for spatially resolved observations. We improve existing models 
and push further our understanding of the disk around HD 1 13766 A. We find that the system is best described by an inner disk located 
within the first AU, well constrained by the Midi data, and an outer disk located between 9-13 AU. In the inner dust belt, our previous 
finding of Fe-rich crystalline olivine grains still holds. We do not observe time variability of the emission features over at least a 8 
years time span, in a environment subjected to strong radiation pressure. 

Conclusions. The time stability of the emission features indicates that /jm-sized dust grains are constantly replenished from the same 
reservoir, with a possible depletion of sub- /jm-sized grains. We suggest that the emission features may arise from multi-composition 
aggregates. We discuss possible scenarios concerning the origin of the warm dust observed around HD 1 13766 A. The compactness of 
the innermost regions as probed by the Midi visibilities, as well as the dust composition, suggest that we are witnessing the outcomes 
of (at least) one collision between partially differentiated bodies, in an environment possibly rendered unstable by terrestrial planetary 
formation. 
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1. Introduction 

Debris disks are byproducts of the planet formation process, 
once the massive, gas-rich circumstellar disk has dissipated. Af- 
ter the gas disks are lost, planetary systems including asteroid 
and comets may remain and evolve through gravitational inter- 
actions. The presence of dust, that originates from collissions 
or evaporation of these bodies, translates into a mid- to far- 
infrared (IR) excess, as demonstrated for the first time for Vega 
(Au mann et al.|1984 i, using the Infrared Astronomical Satellite. 



Zodiacal dust - circumstellar matter - Infrared: stars - Techniques: 



Departure from photospheric emission at these wavelengths is 
the consequence of thermal emission arising from dust grains 
heated by the stellar radiation. Several hundreds of stars are now 
known to harbor such excess emission (e.g., Chen et al.|2 006). 
Our understanding of their nature has progressed during the past 
decades, both observationally and theoretically (see Wyat t|2008| 
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* Based on Herschel observations, OBSIDs: 1342227026, 
1342227027, 1342237934, and 1342237935. Herschel is an ESA 
space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA. Based on Visir observations collected at the VLT (European 
Southern Observatory, Paranal, Chile) with program 089.C-0322(A) 



and |Krivov||20 1 0| for recent reviews). In most cases, these sys- 
tems compare very well with the Kuiper belt observed in our 
Solar System (e.g ., |Carpenter et~aLl[2U09l |Liseau et alT1[2t3T0] 
|Lohne et al.|2012| i, where excess emission is detected in the far- 
IR. Several studies made use of high angular resolution instru- 
ments to detect hot exozodiacal dust around well-known debris 
disks (e.g., Fomalhaut and Vega 
2011 respectively) 
y small (~ 1 



etal. 



ative 



see 

In such cases 



Absil et al. 2009, Defrere 



the near-IR excess is rel- 
h) and departure from photospheric emission 



only becomes significant at mid- or far-IR wavelengths. A small 
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sample of the debris disk population is referred to as "warm de- 
bris disks" (e.g., Moor et al.|2009| > where emission in excess is 
significant in the mid-IR, suggesting the presence of warm dust 
grains, close to the star. Within this class of warm debris disks, 
a few objects are of particular interest since they also display 
emission features in mid-IR spectroscopic observations. These 
emission features are associated to (sub-) /mi-sized dust grains, 
that are warm enough to contribute significantly in the mid-IR 
(see Henning|20 1 0| for a review on optical properties of silicate 
dust). 

These warm debris disks with emission features around 
Solar-type stars are rare (~ 2% of all debris disk s, |Bryden et aL] 
2006] |Chen et al.|2006) . As discussed in |Wyatt| ( |2008| l, the ab- 
normalous brightness of these debris disks cannot be solely ex- 
plained by a steady state collisional evolution. If it was, the colli- 
sional lifetime of the parent bodies responsible for the dust pro- 
duction, and hence the infrared luminosity, would be too short 
compared to the age of the host star. Consequently, two sce- 
narios are possible, either (z) a recent catastrophic collision be- 
tween two planetesimals, or (z'z) an outer belt feeding the inner 
regions, for instance via scattering of bodies by one ore more 
planet ( Bonsor et al.||2012" I or in a similar scenario to the Late 
Heavy Bombardment phase, that happened in the Solar System 
( Gomes et al.||2005 i. The overall scarcity of the "warm debris 
disk" phenomenon would suggest a relatively low probability for 
such events to happen and produce enough dust to be detected. 
Thorough studies of such objects can therefore improve our un- 
derstanding of the final steps of planetary formation, as well as 
provide key clues about the origin of zodiacal dust or Kuiper-belt 
dust in our solar system. 

HD 1 13766 is one of these rare warm debris disks that show 
prominent emission features in the mid-IR. The system is a bi- 
nary of two almost identical stars with spectral types F3/F5, 
with a projected separation of 1.3 " (~ 160 AU at d* = 123 pc, 



Table 1. Final products of Herschel/PACS data processing. "HP" corre- 
sponds to the HighPass filter widths in units of number of readouts. 



van Leeuwen|20 07 ), the dust being located around the primary 



HD 113766 A (L+ = 4.4L Q , psseetal 



spectrum of the 10-16 Myr old debris disk was first studied by 
Lisse et al. ( 2008| l. In Olofsson et al. ( 20 1 2 1 we presented a 



2008). The Spitzer/lRS 



first model to reproduce the Spitzer/lRS spectrum of this source. 
In the following we present new Herschel/PACS photometric ob- 
servations at 70, 100 and 160 /mi, supplemented by archival 
VLTI/Midi interferometric measurements as well as (new and 
archival) VLT/Visir observations in the N- and 2 _rj ands. In 
section|2]we describe the observations and data processing, in 
Sect. [3] we present our model and results. We discuss the impli- 
cations of our findings in Sect.|4]before concluding. 



2. Observations and data processing 

2.1. Herschel/PAcs observations 



We observed HD 1 13766 A with the Pacs photometer (Poglitsch 
|et al |2010] > onboard Herschel ( [Piibratt et al.|2010| >. The dataset 
consists of 4 Pacs mini-scan maps (program "OT1 _jolofsso_l"), 
two in the blue filter (70 /mi, OBSIDs 1342227026 and 
1342227027), and two in the green filter (100 /mi, OBSIDs 
1342237934 and 1342237935). The observing duration was 
220 x 2 sec in the blue filter and 670 x 2 sec in the green fil- 
ter. When observing with Pacs with the blue camera (blue or 
green filters, at 70 and 100 /mi, respectively), data is also simul- 
taneously acquired with the red camera (160 /mi). We therefore 
combined the 4 mini-scan maps in the red filter to obtain the 
photometry at 160 /mi. 



Filter 


Wavelength 


Pixel size 


HP 


Flux 




M 


["] 




[mJy] 








15 


386.5 + 7.2 


ill.. 
Blue 


70 


1.1 


20 


390.1 + 8.7 








35 


391.4 + 9.7 








20 


223.4 + 5.6 


Green 


100 


1.4 


30 


225.4 + 6.6 








45 


226.7 + 7.5 








35 


94.3 + 11.1 


Red 


160 


2.1 


45 


99.3 + 10.1 








70 


113.6+12.6 



The data were processed using the Herschel interactive pro- 
cessing environment (HIPE, version 10.0.667), using standard 
scripts that include bad pixel flagging, detection and correction 
for glitches and HighPass filtering in order to remove the 1// 
noise (see Poglitsch et al.|2010] for more details). We discarded 
the first column of each detector sub-array in order to eliminate 
electronic crosstalk. Finally, we processed the data with different 
HighPass filter widths in order to test the robustness of the data 
reduction. The different values are tabulated in Table[T] (column 
"HP"). 

The aperture photometry in all three filters was extracted 
within HIPE from the combined maps. The aperture radii are 
set at 10" for the blue and green filters and 20" for the red filter. 
The final flux values do not depend on the aperture sizes since 
an aperture correction is performed to compensate for flux loss 
outside of the aperture. This correction is based on a calibrated 
encircled energy fraction parametrisation contained in the cali- 
bration data that is included in HIPE. In order to estimate the 
noise level and thus the quality of the detection in all filters, we 
randomly placed apertures within a distance of 100, 75, and 50 
pixels (for blue, green, and red filters, respectively) to the cen- 
tral pixel. Doing so, we avoided borders of maps that are more 
noisy due to a smaller detector coverage, and are therefore not 
representative of the noise close to the source. We extracted the 
photometry for all of these apertures and then fitted a Gaussian 
profile to the flux distribution histogram to estimate the <x un- 
certainty from the width of the Gaussian. Overall, these uncer- 
tainties are consistent with predictions from the HerschelSpotQ 
software. Figure[T]show the Pacs images at 70, 100, and 160 /mi 
(left to right, respectively), where the source is clearly detected 
at all three wavelengths. Table[T] summarizes the final results of 
the data processing, with pixel sizes for the three filters, aper- 
ture corrected fluxes for the different HighPass filter widths, and 
estimated 1 <x uncertainties in units of mJy. The fluxes used in 
the analysis are the averages of those reported in Table[T] For all 
three filters, the photometric values are almost identical for the 
different HighPass filter widths, within the uncertainties. One 
should notice that the measured 70 /mi flux is consistent with the 
Mips observations reported by |Chen et al.|pOTT[ Fjo =390 mJy 
± 1 1 mJy). 



2.2. Complementary observations 

In addition to the far-IR observations, we included the 
Spitzer/lRS spectrum of HD 1 13766 A (between 5 and 35 //m) in 



http://herschel.esac.esa.int/Tools.shtml version 6.2.0 
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Fig. 1. Herschel/PACS observations at 70, 100, and 160pm (left to right, respectively). 



the analysis (see Ol ofsson et al.|2012 for details on the data pro- 
cessing). We supplemented these data with VLT/Visir observa- 
tions (our program and data published in |Smith et al. 2012|l and 
VLTI/Midi interferometric measurements (S mith et al.|2012| l. 

On the night of the 2012-Apr-21, HD 11 3766 A was ob- 
served with the VLT/Visir instrument (program 089.C-0322(A)) 
in the /V-band. We defined four spectral windows, centered at 
8.8pm (8.0 < A < 9.6pm), 9.8pm (9.0 < A < 10.6pm), 11.4pm 
(10.4 < A < 12.5 pm), and 12.4pm(11.4 < A < 13.5 pm), in low 
resolution mode (R ~ 350 at 10pm). We additionally retrieved 
the spectroscopic observations presented and analysed in Smith 
et al.| ( |2012[ ). These observations were performed on the night 
of the 2009-May- 1 1 in low resolution mode as well, for two 
spectral windows centered at 8.8 and 11.4pm (program 083. C- 
0775(E)). Both datasets were processed using standard scripts 
provided within the ESO Gasgano softwar^] Right panel of Fig - 
ure|2]shows the Irs and Visir spectra obtained in 2004, 2009, and 
2012 (open black circles, thick red, and thin blue lines, respec- 
tively). One should note that both Visir spectra, taken almost 
three years apart agree surprisingly well. However, the spectrum 
obtained in our program in the 9.8 pm spectral window displays 
an offset of about 10% higher in flux. Such "stiching" problems 
are not unusual, but the exact causes are not clear. We checked 
the atmospheric conditions during our observing program: the 
seeing and atmospheric pressures appear to be constant over the 
run. Besides the 9.8 pm band offset, both Visir spectra show a 
deficit of flux of about 30% compared to the Irs spectrum. Smith 
et al. (2012 1 already discussed in detail these differences (e.g., 
spatially extended emission) and concluded these differences to 
be within the absolute calibration uncertainites of Visir (see dis- 
cussion in Geers et al. 2007}. 

The archival Midi data presented in |Smifh et al.| (|2012) were 
processed using scripts dedicated to observations of faint targets 
with Midi (see Burtscher et a l.|2012| for more details) interfaced 
with the MIA+EWS softwar^f] Tabled summarizes the observ- 
ing log for the interferometric measurements. As in Smith et al.| 
(2012), we found that the observations performed on the night 



of the 2007-May-30 (program 079.C-0259(F)) had a low signal 
to noise ratio, and we did not include them in our analysis. The 
calibrated visibilities and correlated fluxes for HD 1 13766 A are 
displayed in the left and middle panels of Fig. [2] The Irs spec- 
trum is also shown for comparison to the correlated fluxes. Be- 



cause of the challenge of calibrating accurately the Visir spec- 
tra, we used the Irs spectrum to compute visibilities from the 
correlated fluxes. The calibrated dataset compares very well 
with the data reduction of Smith et al. ( 2012| l, except for one 
single run (083.C-0775(B), observed at 05:32:55). When cal- 
ibrating these observations with their corresponding calibrator 
(HD 112213, observed at 05:15:45), we found a correlated flux 
much higher than the Irs spectrum. We constructed the transfer 
function based on all the calibrators observed during that night. 
At the time of the science observations, we found a variation 
of about a factor 2 compared to averaged transfer function over 
the rest of the night. This suggests that the observing condi- 
tions were strongly varying. We reduced these science observa- 
tions with a later observation of the same calibrator (06:34:00). 
One can see in Fig. [2] that both visibilities and correlated fluxes 
are large for this run (orange dashed line) and that there is a 
fairly significant difference with the consecutive observations on 
the same night, with the same baseline (06:09:53, orange dot- 
ted line). As explained in Burtscher et al. ( 2012[ ) uncertainties 
are given by the quadratic sum of the photon noise error de- 
livered by the MIA+EWS software and the standard deviation of 
the transfer function over the night of observations. Overall, the 
final uncertainties are dominated by the latter one (between 5 
and 10%). For the "problematic" observing run (083.C-0775(B), 
time 05:32:55), the standard deviation of the transfer function 
over the night is closer to 25%. 



2 http://www.eso.org/sci/software/gasgano/ version 2.4.0 

3 http://home.strw.leidenuniv.nl/~jaffe/ews/ version 2.0 



We did not reduce the VLT/Visir images analyzed in Smith 
|et al.| ( |2012[ ). We instead used the full width at half maximum 
(FWHM) presented in their study. In none of the N- and £>-band 
images (A c = 11.85 and 18.72pm with AA = 2.34 and 0.88 pm, 
respectively) presented by the authors, the debris disk was spa- 
tially resolved compared to a reference PSF They fitted the az- 
imuthally averaged surface brightness distributions with Gaus- 
sian profiles with FWHM of 0.322 " and 0.498 " in the N- and 
2-bands, respectively. During the analysis, we will check if syn- 
thetic images convolved with such PSFs are consistent with the 
measurements of |Smith et al.| ( [20T2] l. After convolving the syn- 
thetic images with the reference PSF, we will compare the sur- 
face brightness distribution to check for possible resolved emis- 
sion in both bands. 
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Fig. 2. VLTI/Midi visibilities and correlated fluxes for HD 1 13766 A (left and middle panels, respectively). The baselines used and their lengths 
are reported in the caption of the left panel. Right panel shows the Visir ,/V-band spectrum at different epochs (2009 and 2012, in thick red and 
thin blue, respectively). The Spitzer/lRS spectrum is shown in open black circles in the middle and right panels. 



Table 2. Log of the VLTI/Midi observations and averaged visibilities in the range 10.5 /urn ± 0.5 /urn, with 1 a uncertainties (see text for details, 
Sect.[3](. The (u, v) values are in m. 



Program 


Date 


Time 


Baseline 


(h,v) 


V 


079.C-0259(G) 


2007-Apr-08 


01:11:11 


UT1- 


-UT3 


(4.9,101.7) 


0.48±0.03 


079.C-0259(G) 


2007-Apr-09 


00:33:40 


UT1- 


-UT3 


(-4.7,101.7) 


0.45+0.03 


083.C-0775(D) 


2009-May-08 


05:15:29 


UT1- 


-UT3 


(65.3,49.9) 


0.61±0.03 


083.C-0775(D) 


2009-May-08 


06:11:31 


UT1- 


-UT3 


(61.8,38.7 


0.68+0.04 


083.C-0775(B) 


2009-May-09 


05:32:55 


UT2- 


-UT4 


(67.2, -33.0) 


0.75+0.10 


083.C-0775(B) 


2009-May-09 


06:09:53 


UT2- 


-UT4 


(57.0, -40.3) 


0.67+0.04 



3. Multi-technics, multi-wavelengths modelling 



As in Olofsson et al. ( 2012| >, a given debris belt is fully de- 



scribed by four free parameters plus a dust model: the inner 
and outer radii of the dust belt (r; n and r out ), a volume density 
exponent a (E(r) = 2o(r/ro) cf , a < 0) and an exponent p for 
the grain size distribution (dn(s) oc s p ds, p < 0). The So refer- 
ence value is found by adjusting the total mass of the dust M^ ust , 
which in turn is obtained by fitting the SED. The relative abun- 
dances for the different dust species are also free parameters in 
the modelling and are found by fitting the SED and the emis- 
sion features detected in the Irs spectrum. For each set of n n , 
r out , a and p, the fitting of the abundances is performed via a 
Levenberg-Marquardt algorithm on a linear combination of the 
thermal contributions of all the dust species considered. The 
best-fit disk model is found via a Monte-Carlo Markov Chain 
(MCMC) exploration of the parameter space (emcee package, 
Foreman-Mackey et al. 2012[ ). All the uncertainties quoted in 
the next sections represent the standard deviations of all the val- 
ues tested during the fitting process. 

In the following, we will first present a model of the debris 
disk assuming that the dust is distributed in a single belt around 
the central star (referred to as model "IB" for "one belt"). We 
will test the robustness of this model, comparable to the one pre- 
sented in Olofsson et al. (2012 1, with respect to the new far-IR 
photometric observations and mid-IR interferometric measure- 
ments. The Herschel data provide new constraints on the cold 
counterpart of the debris disk, while the mid-IR Midi data are 
more sensitive to the warmer dust and provide direct information 
on the geometry of the dust distribution. As a matter of fact, the 
interferometric observations suggest the region emitting in the 
iV-band is compact (< 10 mas, or 1.2 AU in diameter, according 
to |Smith et al.||2012[). On the other side , our previous best-fit 
model described in Olofsson et al.|2012 was found to be fairly 
extended (r out ~ 12 AU, with a density profile with an exponent 



of a = -1.8). We can therefore already anticipate that fitting the 
SED and the Midi data simultaneously will be challenging. 

In order to compare our model to the Midi and Visir obser- 
vations we updated the Debra code to produce synthetic images 
at different wavelengths. To do so, at a given wavelength, we 
distribute the thermal surface brightness distribution (F v (r), be- 
tween r m and r out ), assuming the disk is axisymmetric and in- 
finitely flat. Possible inclinations effects are taken into account 
by projecting pixels of the image onto the sky plane. However, 
because we assumed the disk to be infinitely flat, we chose the 
inclination i (0° being face-on) to be smaller than 60°. In the 
case of HD 1 13766 A, this choice is supported by the analysis of 
Smit h et aI7| ( |2012j > who could successfully model the Midi data 
with a circularly symmetric flux distribution (i.e., small i). As 
we are solely interested in the variation of the visibilities as a 
function of the baseline lengths, we averaged the Midi visibil- 
ities in the spectral range 10.5+0.5/im for each Midi baseline. 
We tested different values for i and the position angle (PA) of 
the debris disk. We tested four different values for i (0, 20, 40, 
and 60°) and 30 values linearly distributed between and 180° 
for the PA (0° meaning the semi-major axis is oriented along the 
North-South direction). 

3.1. One belt (model "IB") 

At first, we tried to model solely the SED, including the new 
Pacs photometric points, and the Irs spectrum. As mentioned 
earlier the Midi data suggest that the regions emitting in the TV- 
band is compact (r < 1.2 AU). Therefore, we first want to test 
how robust the assumption of a single dust belt is, with respect 
to the far-IR measurements, and underline at the same time the 
limitations of such approach. We will then compare the best-fit 
"IB" model to the interferometric data and discuss the subse- 
quent implications on the spatial distribution of the dust. The 
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Fig. 3. Upper and lower panels: Results for the "IB" and "2B" models. Left panel: SED of HD 1 13766 A. The Irs spectrum is in black, Pacs 
data are shown as black open circles. For the "IB" model, the final best fit is in red. For the "2B" model, the final total fit is in green, the separate 
contributions of the inner and outer belts are in red and blue, respectively. Middle left panel: blow-up of the stellar subtracted Irs spectrum. Color 
coding is the same as for the left panel. Residuals (100 x [F v obs - F v raod ]/F Vi0bs ) are displayed below the spectra. Middle right panel: Observed 
and modelled Midi visibilities as a function of the baseline (open black circles and red squares, respectively, with 1 cr uncertainties). Right panel: 
modelled surface brightness at 10 and 20yum (thin and thick black lines, respectively) of synthetic images convoled with 2D Gaussian with FWHM 
of 0.322 and 0.498". The 2D Gaussian with the aforementioned FWHM, assumed to be PSF representative of the images, are shown in thin and 
thick dashed red lines, for comparison. 



dust properties (absorption efficiences and grains sizes) are the 
same as the ones described in Olo fsson et al.| ( |20T2] l. Grains sizes 
for the crystalline dust grains are 0.1 < s < 1 fim and 0.1 //m 
< s < 1 mm for the amorphous dust grains. The minimum and 
maximum values for the grain sizes are not free parameters in the 
modelling. This choice mostly originates from the experimental 
setup for the laboratory experiments that led to the extinction ef- 
ficiencies of crystalline olivine grains (see Olofsson et al.|2012| 
for more details). 

Modelling only the SED with the new far-IR observations, 
we now find the dust to be located between (0.4 + 0.2) AU 
and (51.9 + 25.2) AU, with a volume density exponent of a = 
-1.9 ± 0.3, and an exponent for the grain size distribution 
p = -3.0 ± 0.1. The total dust mass is 1.3 x 10~ 2 M ffi . The 
relative abundances for the different dust species are as fol- 
lows: (43.1 ± 1.4)% of amorphous silicates (with a predom- 
inance of MgFeSi206 and MgFeSi04 grains, optical constant 
from |Dorschner et al.|1995] ), (17.3+0.7)% of Mg-rich crystalline 
olivine grains, (2 1 .0 + 0.6)% of Fe-rich crystalline olivine grains 
(laboratory measurements from Tamanai & Mutschke 201 0|fo r 
the olivine data), (11.6 ± 0.5 )% of enstatite ( [Jager et al.|1998a| l, 
(0.0+0.8)% of/3-cristobalite ( [Tamanai|20f0l>, and o f (7.0+0.6)% 
amorphous carbon dust grains ( |Jager et al.|1998b| mostly in the 
form of cellulose pyrolized at 800°C, "cel800")- Disk parame- 
ters for the best-fit model are reported in the first line of Table[3] 

Fits to the SED, the stellar subtracted Irs spectrum, Midi 
data, and a comparison to the Visir images are displayed on the 
upper four panels of Figure|3](left to right, respectively). For the 
SED, the Irs spectrum is shown in black, the Pacs data are shown 
as black open circles, and the best fit is represented as a solid red 
line. The color coding is the same for the blow up of the fit to the 
Irs spectrum (second column), and residuals are shown below 
the fit. Residuals are computed as 100 x [F Vi0 b s - ^ v ,mod]/^v,obs, 
where F vo b s and F^ moi i are the observed and modelled fluxes, 



respectively. The third column panel shows the observed and 
modelled Midi visibilities for the best-fit to the SED (uncertain- 
ties are 1 cr). Right panel shows the comparison between the 
PSFs fitted to the Visir observations ( [Smith et al.|2012| l, in the 
TV- and 2-bands and the azimuthally averaged profiles of syn- 
thetic images for the final best-fit model, convolved with PSFs 



with FWHMs of 0.322" and 0.498", respectively, as in |Smith 
|et^Ll ( |20T2l l. 

While the best-fit model agrees well with the SED and the 
Irs spectrum, the observed Midi visibilities are severely under- 
predicted. As discussed above, this can be explained by the rela- 
tive compactness of the emitting regions probed by the interfero- 
metric data. The A^-band emission in our best-fit model is too ex- 
tended and therefore the modelled visibilities drop very quickly. 
This is also true for the Visir observations. The modelled sur- 
face brightness profiles over-predict the observed 2D Gaussian 
profiles assumed to be representative of the PSFs, in the 7A/-band 
and even more in the Q-band (the disk was not resolved in the 
Visir observations). 

As a final test, we tried to include the Midi observations in 
the fitting process (by summing the reduced xl from the fits to 
the SED and Midi data), but could not achieve any better results. 
Both the Midi and the mid- to far-IR part of the SED are under- 
predicted. This can be explained as a tradeoff between the com- 
pactness of the dust belt and its contribution in the far-IR. Under 
the assumption of a single dust belt, the fits to the SED and in- 
terferometric data are diametrically opposed. From a modelling 
point of view, the main conclusion of this exercise is that SED 
modelling can be misleading: as shown in the upper left panel of 
Fig. [3 we could obtain a decent fit to the SED, but this model is 
quick y proven wrong as soon as spatially resolved observations 
are available. 
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Fig. 4. Probability distributions from the MCMC run for model "2B", 
as a function of relative abundances of Mg- and Fe-rich olivine grains 
(solid black, and dashed red, respectively). 



3.2. Two belts (model "2B") 



Several recent studies have shown that debris disks may contain 
several planetesimal belts (e.g., Morales et al. 2011 or the disk 
around r Ceti, |di Folco etaLl|2007[ |Greaves et al.||2004| >. We 
have shown in Sect. 3.1 that the assumption of one single dust |et al. ( |2012 



at larger distances (Sect. |4.1] i. We set the minimum grain size to 
■Smin = 4 pm for all dust species in the outer belt, as there are 
no observational evidences, such as emission features at long 
wavelengths, that would require including smaller dust grains. 
Because we assumed in our model that crystalline dust grains 
have sizes smaller than 1 pm, they are not accounted for in the 
outer belt (see Sect. |4.3| for further discussion). This time, for 
each model of the MCMC run we fitted the SED and the Midi 
observations simultaneously. 

For the inner dust belt, we obtained a best fit to the data with 
dust being located between (0.6 + 0.1) AU and (0.9 ± 0.2) AU, 
a density profile in a = —1.7 ± 0.3, and a grain size distri- 
bution exponent of -3.5 ±0.1 (consistent with the prediction 
for a collisional cascade p — -3.5, Dohnanyi 1969). The to- 
tal dust mass is 9.0 x 10~ 5 M @ . The relative abundances for 
the inner belt are the following: (44.5 ± 1.4)% of amorphous 
silicate grains (dominated by Mg2SiC<4 grains, optical constant 
from Jage r et al.|2003] >, (17.0 ± 0.8)% and (18.8 + 0.8)% of Mg- 
rich and Fe-rich crystalline olivine grains, respectively, as well 
as (7.7 ± 1.2)% of enstatite, (9.4 ± 0.4)% of /3-cristobalite, and 
(2.7 + 0.5)% of amorphous carbon grains (cellulose pyrolized at 
800° C, "cel800"). Figure|4] shows the marginalized probability 
distributions built from the MCMC run, for the relative abun- 
dances of Mg- and Fe-rich crystalline olivine grains (solid black 
and dashed red histograms, respectively). Even though more Fe- 
rich crystalline olivine grains are found compared to Mg-rich 
olivine, both distributions are overall close to each other. This 
result differs from the abundances previously found in Olofsson 



belt around HD 113766 A cannot account for all the available 
dataset. We can therefore start discussing the possibility of two, 
spatially separated dust belts around the central star: an inner 
dust belt, compact, close to the star that can reproduce the Midi 
data and the near- to mid-IR flux, and a second dust belt, far- 
ther away from the star that will contribute to the emission in 
the far-IR. To account for an additional belt, the modification to 
be applied to the Debra code is straightforward as we simply 
duplicated the computations done for one dust belt to the sec- 
ond one. A second radial grid, with different r; n and r out values, 
with its own volume density (a) is generated, as well as the cor- 
responding dust temperatures as the function of the grain sizes 
and compositions. The dust properties in the outer belt can be 
defined separately from the ones in the inner belt (different dust 
species, s m [ n , s max , and p). Finally, the fitting to the SED is per- 
formed over the sum of the contributions of the two dust belts 
(returning the dust mass for each belt). 

As discussed in Olofsson et al.| ( |2012"| ), the radiation pressure 
force is strong around HD 1 13766 A. We computed the unitless 
j6 ratio between the radiation pressure and gravitational forces, 
for different grain sizes and dust compositions as described in 
~( [T979) . 



Burns et al. 



We find that grains with sizes smaller than 
(4 + 0.5) pm (depending on their compositions) are expected to 
be short-lived around the central star (see Sect. |4. 1} . However, 
the detection of strong emission features indicates that warm 
(sub-) yum-sized grains must be located close to the star to con- 
tribute in the Irs spectrum. Therefore, the dust properties (ab- 
sorption efficiencies and grain sizes) of the innermost dust belt 
in model "2B" (as in "two belts") are chosen to be similar to the 
"IB" model (Sect. |3.1| , regardless of the radiation pressure is- 
sue. The relative abundances are still free parameters. For the 
outer dust belt, we did not include such small grains. The ra- 
tio ft between the radiation pressure and gravitational forces is 
distance-independent (both forces decreasing as r~ 2 , their ratio 
no longer depends on r), while the collisional timescale is longer 



4 times more Fe-rich olivine grains than Mg-rich 
crystalline grains). This difference can be explained as the outer 
planetesimal belt acts as a source of featureless, continuum emis- 
sion. Consequently, this modifies the required contribution of the 
inner dust belt that is responsible for the emission features, and 
hence the relative abundances of the dust species. Nonetheless, 
it means that Fe-rich crystalline olivine dust grains are neces- 
sary to provide a good fit to the emission features detected in the 
mid-IR spectrum. 

For the outer belt, we found r m = (9.0 + 1.4) AU and r out = 
(13.1 ± 3.4) AU, with a = -1.2 ± 0.3, and p = -3.8 ± 0.2. The 
dust mass in the outer planetesimal belt is of about 2.4xl0~ 3 M m . 
The dust content of the outer planetesimal belt consists entirely 
of amorphous carbon grains (68.7% of "cel600" and 31.3% of 
"cellOOO"). This can be explained by our choice of s m [ n - 4pm 
which damps most of the emission features, as well as the lack 
of significant emission features in the spectral range where the 
outer dust belt contributes most. Similar results were found 
by |Lebreton et al.| ( |2012[ ) in the process of modelling the SED 
of HD 181327. When considering only silicates and carbona- 
ceous grains, carbon grains were preferentially used to fit the 
data. The authors pushed their models further by including the 
effects of porosity as well as the contribution of icy grains, re- 
ducing at the same time the relative fraction of carbon grains. 
For HD 113766 A, given the intertwined contributions of both 
the inner and outer belt, and the limited spectral coverage at 
wavelengths longer than 160 pm, we did not perform similar in- 
vestigations. One should also note that our result for the location 
of the outer belt is highly dependent on the minimum grain size 
■Smin opted in our model. We find typical temperatures for the 
carbon grains between 100 and 200 K. To have comparable tem- 
peratures, if s m j n > 4 pm, grains will have to be located closer to 
the star, and vice-versa if s m j n < 4 pm. Spatially resolved obser- 
vations of the outer belt are mandatory to better constrain its lo- 
cation and the dust content. Even though the disk is not detected 
in the Visir images, they provide an upper limit of about 0. 14" 
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Table 3. Best-fit results for the disk parameters for both the "IB" and "2B" models 



Model 


Belt# 


fin 

[AU] 


r out 

[AU] 


a 


P 


M dust 
[M e ] 


One belt ("IB") 


1 


0.4 ± 0.2 


51.9 ±25.2 


-1.9 ±0.3 


-3.0 ±0.1 


1.3 x 10~ 2 


Two belts ("2B") 


1 

2 


0.6 + 0.1 
9.0 ± 1.4 


0.9 ± 0.3 
13.1 ±3.4 


-1.7 ±0.3 
-1.2 ±0.3 


-3.5 ±0.1 
-3.8 ±0.2 


9.0 x 10" 5 
2.4 x 10" 3 



(17 AU at 1 23 pc) for the maximal extent of the disk ( Smith et al. 
|2012) ). Disk parameters for the best-fit model are reported in the 
last two lines of Table|3] 

The lower four panels of Fig.[3]show the results for the "2B" 
model, with the same arrangement as described in Sect. |3.1| For 
the left and middle left panels (SED and Irs spectrum) the to- 
tal fit is shown as a solid green line, while the two contribu- 
tions of the inner and outer dust belts are represented as red and 
blue solid lines, respectively. Compared to the "IB" model, the 
best-fit model is much more compact and therefore both the Midi 
and Visir data are much better reproduced (except the "problem- 
atic" Midi run discussed in Sect.[2]i. As shown by the SED, the 
contribution of the outer belt in the A^-band is almost negligi- 
ble compared to the total flux. Therefore, the Midi observations 
are mostly sensitive to the inner belt, which is much more com- 
pact. We find its inclination to be of 40° with a position angle of 
120°, however, one should note that the inclination is sampled by 
steps of 20° only. The overall compactness of the inner dust belt 
also explains why the modelled A^-band Visir synthetic image 
agrees with the observations, as most of the emission at 10/itri 
arises from the inner belt. As for the Q-band image, [Smith et aL 



pOT2| l found an upper limit at 0.14" (17 AU at 123 pc) for the 
maximum extent of the disk, consistent with our r out =13 AU. 

Our best fit model indicates that the inner and outer belt has 
to be spatially separated. However, we further pushed our mod- 
eling effort and tried to fit the dataset with a single dust belt 
that shows a "break" in the density distribution (i.e., a single 
but not continuous belt as opposed to Sect ]3. 1[ >. The motivation 
being to mimic an inward radial drift of dust grains, followed 
by the formation of a dust ring at the sublimation radius, as de- 
scribed in Kobayashi et al. ( |2009 1 for instance. We therefore 
run a new set of models, for which we force the inner radius 
of the outer belt to be the same as the outer radius of the in- 
ner belt (seven free parameters instead of eight). The best fit 
results for the SED and the Midi data could not reach the qual- 
ity of the best-fit "2B" model. We find that the disk extends 
from 0.6 AU up to 16.5 AU, with a break in the density at about 
0.7 AU. The fit to the Irs spectrum is acceptable except for the 
red wing of the 10/im emission feature where the modeled flux 
is larger than the observed flux. Because the disk extends far- 
ther out than ~ 1 AU, the interferometric measurements are still 
under-predicted (V ~ 0.5 - 0.6). To compare our results with 
the calculations presented in Kobayashi et al. ( 2009 ), we com- 
pute the number density n as a function of the distance r. For 
the best-fit model, the break at the intersection between the two 
adjacents dust belts has to be extremely large: the quantity n x r 
is more than five orders of magnitude larger in the innermost re- 
gions (r < 0.7) compared to the outermost regions of the disk 
(r > 0.7). If the warm ring was the consequence of radial drift 
and dust sublimation, Kobayashi et al. (2009) predicted an en- 
hancement factor smaller than 10 for silicate dust grains around 
a F-type star. We can therefore confidently rule this scenario out 
to explain the origin of the transient dust around HD 1 13766 A. 



As mentioned in Sect.[T[ the secondary of the binary system 
is located at 1.3" of HD 113766 A (160 AU at 123 pc). Without 
a fully characterized orbital solution for the binary system it is 
challenging to discuss in detail the influence of the secondary on 
the debris disk. However, it is unlikely that it can have a ma- 
jor effect on the outer disc located at ~ 13 AU from the primary. 
The only possibility would be if the companion star was on a 
very eccentric orbit with a periastron passing close to the outer 
belt, but in this case the secondary would spend most of its time 
in remote regions close to its apoastron, leaving much time for 
the circumprimary disc to "recover" from the perturbations at 
periastron passages (Thebau lt et al.|2010 Thebault|2012| l. 



4. Discussion 

Results presented in Sect.[3] show that the available dataset for 
HD 113766 A can be best reproduced under the assumption of 
two spatially separated dust belts. Even though the SED alone 
can be reproduced under the assumption of a single dust belt, 
the Midi observations proved extremely useful to constrain the 
spatial extent of the region responsible for the emission in the 
A^-band. If the inner dust belt extends farther out than ~ 1 AU, 
the modelled Midi visibilities are greatly under-predicted at all 
baselines. Our model is in good agreement with the results of 
Smit h et"aT] ([2012) who modelled the Midi data with either a 
symmetric Gaussian model with a FWHM of 10 mas (1.2 AU 
in diameter) or a ring model with a radius of 6 mas (~ 0.7 AU). 
Our modelling results, supported by the new Pacs observations 
enabled us to go one step further and demonstrate the presence 
of an additional planetesimal belt farther out in the system. 

4.1. Time variability 

Even though the absolute calibration of the Visir spectra is chal- 
lenging, one can already see from the right panel of Fig|2]that the 
shape of the emission features has not changed between 2004, 
2009, and 2012. The relative strengths, the widths, and peak po- 
sitions of the features seem to remain the same over a time span 
of at least 8 years. Such a phenomenon appears to be also valid 
for other debris disks with strong emission features (Olofsson 
et al.|2012 Johnson et al.|2012 1, while rapid dissipation of cir- 
cumstellar material has only been observed for one single object 
(Meli s et al.|2012| l. If (sub-) yum-sized crystalline dust grains are 
solely responsible for the emission features observed at 10 and 
1 1 //m, it would suggest that this population of grains is replen- 
ished from the very same reservoir (same composition), on the 
timescale of the evacuation by radiation pressure. To compute 
this timescale f, p , we followed the formalism to be presented in 
Lebreton et al. (in prep, private communication). The authors 
define f lp , which depends on the ft ratio, as the time needed for 
a dust grain to be transported from its initial distance r to twice 
this distance 2xr. As mentioned earlier, we computed the ft ratio 
for grains of different compositions as a function of their sizes. 
For HD 1 13766 A, this ratio peaks at ft — 3 for grains with sizes 
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of about 0.5 /mi. For these dust grains, assuming r = 0.75 AU 
(l/out - nn]/2), we find f rp ~ 1 month. On the other hand, small 
dust grains are continuously produced by collisions. To estimate 
the collisional lifetime of the smallest grains in the inner dust 
belt, we used the formula f c = (Qt)~\ where £2 is the angular 
velocity at a distance r and r is the geometrical vertical opti- 
cal depth (Thebault & Augereau 2007). Even though the exact 
collisional timescales may depart from this value, it provides an 
order of magnitude estimate. At a distance r = 0.75 AU, we find 
a collisional lifetime - source of dust replenishment - of about 2 
months, twice the timescale of evacuation by radiation pressure 
frp. For the outer belt, assuming r = 10 AU, we find f rp ~ 5 years 
(for grains with j3 — 3) and f c = 300 years, which supports our 
choice of s m in = 4 /mi in the outer belt. One should note that the 
timescale for the Poynting-Robertson drag, is much longer com- 
pared to the aforementioned timescales and can be neglected in 
the case of HP 113766 A (Wyatt|2005| l. 

The sub- //m-sized grains responsible for the observed emis- 
sion features should be evacuated rapidly from the inner belt. 
One may therefore consider dust aggregates of different com- 
positions (e.g., Lindsay et al. |2010[ ) as a possible explana- 
tion. Computations of absorption coefficients, making use of 
the discrete-dipole approximation method (Draine 1988) have 
shown that small crystalline monomeres may be incorporated 
into larger dust aggregates but still show prominent emission fea- 
tures associated with crystalline dust grains (e.g., |Wooden et aE] 
|2012| l. These larger aggregates would then be less subjected to 
radiation pressure and stay gravitationally bound to the central 
star on longer timescales. But even these aggregates will even- 
tually collide, fragment, and be blown out of the system. This 
suggests that replenishment from a reservoir with the same bulk 
composition has to occur. Another way to explain the overall sta- 
bility of the emission features would be to introduce a deficit in 
the grain size distribution for grains that are affected by radiation 
pressure (e.g., |Johnson et al.|2012| l. However, this modelling is 
out of the scope of this paper given the total number of free pa- 
rameters and the lack of spatially resolved observations to better 
constrain the location of the dust. 



4.2. The optical depth and opening angle 



Natively, the model space in the Debra code is ID: there is only 
a radial dependence on the density and temperatures of the dust 
grains. However, an infinitely flat disk would quickly become 
radially optically thick. Therefore, the vertical dimension is de- 
scribed via an opening angle 9 (0° being an infinitely flat disk). 
At constant volume density, increasing 9 decreases the local den- 
sity in all of the radial cells, and hence the optical depth becomes 
smaller as well. For a given disk model, the Debra code returns 
a threshold value 9\{ m , which is the smallest opening angle for 
which the debris disk is still in the optically thin regime. One 
should note that for small opening angles, 9 is comparable to 
twice the aspect ratio of the disk 2H/r (same formalism as in 
Thebault 2009 ). For the best-fit "2B" model, we found a £>ii m , inner 
of about 3.8° for the innermost belt, a value consistent with the 



4.3. The origin of the transient dust 



one found by Lisse et al. (2008 4°). The opening angle for the 
outer belt is smaller (0 um> outer ~ 0.9 °). The lower limit for the 
opening angle 0ii m> i nne r is in good agreement with the work pre- 
sented in Thebault (2009 ), who found a minimum "natural" ob- 
served aspect ratio of ~ 2.3 ° (defined as 9/2). 



In |Qlofsson et al.| ( |2012| > we hypothesized that the transient dust 
grains observed around HD 1 13766 A may be the outcome of a 
disruptive collision of partially differentiated bodies. This con- 
clusion was based on specific mineralogic markers: the com- 
position of crystalline olivine grains containing a Fe/[Fe + Mg] 
fraction of about 20%. Based on the new observations available, 
we refined our previous best-fit model by including an additional 
planetesimal belt farther away from the star. How does this new 
model fit in this picture ? One scenario that has been suggested 
in the literature to explain the excess emission observed around 
warm debris disks is that an outer belt is "feeding" the inner- 
most regions. Comet-like bodies from the outer regions may be 
scattered inwards, for instance by a (chain of) planet(s) as de- 
scribed in Bonsor et al. (2012). These bodies will then collide 
with each other or sublimate and release small dust grains in the 
immediate vicinity of the star. With our finding of an outer dust 
belt, one may wonder if such an event is taking place around 
HD 11 3766 A. 

The TY-band Midi observations ruled out any significant 
emission to be outside of ~ 1 AU, suggesting that the phe- 
nomenon responsible for the transient dust has to be "local", as 
opposed to extended. For icy comets scattered inwards from the 
outer belt, sublimation may release a population of dusty grains. 
The majority of these grains would be released inside of 3.3 AU, 
based on heating of large, km-sized comets by a stellar radiation 
field from a 4.4L Q star (private communication with Amy Bon- 
sor and Ulysse Marboeuf, Bonsor et al. in prep). The fate of 
such grains released from a comet is unclear- whilst they could 
supply the observed excess emission, release of dust grains from 
an eccentric comet does not by itself explain the compact nature 
of the inner dust belt. If the observed warm material originated 
from the outer regions, the Midi data could only be explained by 
the capture of a massive comet that is out-gasing and evaporat- 
ing close to the star (a local production of dust, Beic hman et al.| 
2"5031 concerning HD 69830). 

Our best-fit model can successfully reproduce the emission 
features detected in the Irs spectrum without including crys- 
talline dust grains in the outer belt. This choice, initially mo- 
tivated by the strength of the radiation pressure can also be justi- 
fied by the lack of emission features associated with crystalline 
grains beyond 30 fim. If the outer belt would contain small crys- 
talline dust grains, they would be cold enough to be respon- 
sible for emission features at long wavelengths (especially the 
33 //m feature arising from crystalline olivine grains). The ther- 
mal emission of the outer belt peaks at about 25-30/mi, and the 
fact that the 33 //m feature is not detected suggests that such //m- 
sized crystalline grains are not present, at a detectable level, in 
the outer belt. The effect of crystalline monomeres inclusions in 
aggregates has been studied for the 10//m emission feature (e.g., 
Min et al. 2008 ), but not for other emission features at longer 
wavelengths. Further investigations have to be performed to es- 
timate how much crystalline grains can potentially be included 
in aggregates that would still not show any detectable emission 
feature around 30 /mi. Overall, this suggests that the crystallinity 
fraction has to be smaller in the outer belt compared to the in- 
ner belt, pointing towards a different composition between the 
two dust belts. Would the mineralogical differences between the 
two belts be confirmed, this would argue against the capture of a 
massive comet originating from the outer regions. 

Our spectral decomposition results suggest that the crys- 
talline olivine grains in the inner belt contains more Fe than what 
is expected for comet-like bodies (Wooden 2008 ; Zolensky et al. 
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2008 ; [Hanner & Zolensky 2010), cold crystalline olivine grains 



observed in other debris disks (J3 Pic, |de Vries et al.||2012 



or 



even younger proto-planetary disks (e.g., Olofss on" et al.||20T0 
Stur m et aLl|2010[ Sturm et al. submitted, and |Juhasz et al.] 



2010 who found an upper limit of 10% for the Fe fraction). 



As previously discussed in Olofsson et al. (20121, crystalliza- 
tion via thermal annealing or gas-phase condensation are un- 
likely to produce Fe-bearing crystalline olivine grains in proto- 
planetary disks ( |Nuth & Johnson|[2006l [Gail|[jo04| . Observ- 
ing a significant fraction of Fe-rich crystals (more than half the 
total abundance of crystalline olivine grains) indicates that the 
dust composition has been altered. Internal metamorphism of 
the olivine content, as discussed in Nakamu ra et al.| ( |201 1) con- 
cerning the asteroid Itokawa, can provide an explanation for our 
finding. Consequently, with present-day observations and labo- 
ratory measurements of optical properties, we reiterate our pre- 
vious conclusion that we may be witnessing the aftermath of 
a disruptive collision within the first AU from the central star. 
The inferred composition suggests that the dust must have been 
processed in a way which is not consistent with our current un- 
derstanding and knowledge about the mineralogy of cometary 
bodies and proto-planetary disks. A collision of differentiated 
bodies as the origin of the transient dust seem to provide a better 
explanation of our results. If so, this conclusion also suggests 
the parent bodies do not originate from the outer belt, support- 
ing the scenario of a "local" event. As discussed in |Triel off & 
Palme ( 2006 ), thermal metamorphism by internal heating is ex- 
pected to be more efficient inside the first planetesimals that are 
formed, as they accrete more short-lived radionuclides ( 26 A1), 
while planetesimals that accreted material later on remain undif- 
ferentiated. As the timescale for planetesimal formation depends 



on the distance to the star (and the surface density, Kenyon & 
|Bromley||2004| |Wyatt||2008| > planetesimals in the outer belt are 
formed later on and thus may remain undifferentiated 



In a recent study, |Jackson & Wyatt (2012 1 investigated the 
outcome of a massive collision in the framework of the forma- 
tion of the Moon. They found that such an impact would re- 
lease enough dust for the corresponding excess emission to be 
detectable in the mid-IR up to 25Myr after the collision hap- 
pened. Their numerical simulations suggest that the orbits of 
the debris released upon the impact would circularize in about 
lOkyr. Naturally, to take place, such event require the presence 
of large planetesimals, implying that km-sized bodies must have 
formed around HD 1 13766 A. As detailled in the review by Wy- 
[att| ( |2008] and references therein), planet formation models sug- 
gest that such bodies had sufficient time to be formed in the 10- 
16Myr old system. In turn, the chaotic phase of planetesimal 
growth in the innermost regions scatters the orbits of km-sized 
bodies and increases the probability of a collision (e.g., [Kenyon] 
& Bromley|2006) . 



5. Conclusion 

In this study, we presented a new model for the debris disk 
around HD 113766 A. This model accounts for new far-IR pho- 
tometric observations as well as archival mid-IR interferometric 
measurements. We underlined the limitations of SED modelling 
and the importance of spatially resolved observations. Our ap- 
proach of modelling simultaneously the SED from near-IR to 
far-IR wavelengths and interferometric measurements has en- 
abled us to have a deeper understanding of the morphology of 
the debris disk around HD 113766 A. We found that the avail- 
able dataset can be best explained by two spatially separated 
dust belts. A first inner dust belt is located within 1 AU from 



the star, a distance well constrained by the Midi observations. 
This inner belt contains warm sub-jum-sized dust grains. These 
small grains are responsible for the emission features observed 
via mid-IR spectroscopy. An outer planetesimal belt, located 
between 9-13 AU (distances much less constrained by the ob- 
servations) account for the far-IR emission observed with Pacs. 
The new observations support our previous conclusion (Olofs- 
son et al.|2012) , where we suggested that the observed transient 
dust is the result of a disruptive collision. This collision could be 
the consequence of on-going terrestrial planet formation around 
the young debris disk (10-16 Myr). 
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